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ABSTRACT 

The jets of GRS 1915+105 carry a considerable energy away from the central source into the ISM. The 
similarity of the jets of this source and jets in radio galaxies or radio-loud quasars suggests that we should 
detect large-scale, synchrotron emitting radio structures surrounding GRS 1915+105. However, these large 
structures have not been found. We show that by adapting a model for the radio lobes of extragalatic jet sources 
we predict a radio surface brightness of the equivalent structures of GRS 1915+105 below the current detection 
limits. The model uses an energy transport rate of the jets averaged over the jet lifetime. This transport rate 
is found to be considerably lower than the power of the jets during the rare major ejection events. Thus the 
lobes contain less energy than would be inferred from these events and produce a lower radio luminosity. The 
model also predicts a lifetime of the jets of order 10 6 years and a gas density of the ISM in the vicinity of GRS 
1915+105 of ~ 150 cm" 3 . The impact sites of the jets are identified with two iras regions with a flat radio 
spectrum located on either side of GRS 1915+105. Observations of molecular lines and dust emission from 
these objects are consistent with our interpretation. Distance estimates for the iras regions give 6.5 kpc and our 
model implies that this is also the distance to GRS 1915+105. This low distance estimate in combination with 
the observed motions of jet ejections on small scales yields a jet velocity of about 0.7 c and an angle of 53° of 
the jets to our line of sight. 

Subject headings: stars: individual (GRS 1915+105), binaries: general, stars: outflows, ISM: jets and outflows, 
ISM: individual (iras 19124+1 106, iras 19132+1035), radio continuum: ISM 



1. INTRODUCTION 

Collimated outflows or jets are by now a commonly ob- 
served feature of X-ray transient binary systems. The sim- 
ilarity of some of the properties of these jets with those of 
the much larger jets in radio-loud active galaxies has led to 
the name microquasars fo r these objects (for a review see 
Mir abel & Rodriguezlll999|) . However, there are also some 
differences between the jets of microquasars and those of ra- 
dio galaxies and radio-loud quasars. One of these is the lack 
of observational evidence for an interaction of the jet flows 
with their environment in the case of most microquasars. 

Discrete ejections of gas from the central binary system 
sweep up the gas in front of them. The associated transfer 
of momentum to this gas should slow the ejections down un- 
less the ejections are much denser than the surrounding gas. 
However, to date only two examples of dece lerating jet ejec- 
tions have been observed for X TE J 1748-2 8 8 ( Hiellming et alJ 
fT99l and for XTE J1550-564 dCorbel et alJl2Q03h . All other 
observed jet ejections fade to below the detection limit before 
any noticeable deceleration has taken place. Although this 
supports the notion of very dense, ballistic ejections in micro- 
quasar jets, such ejecta must also slow down eventually due to 
the continued loss of momentum to the swept up gas in front 
of them. 

Another possibility is that the continued ejection of jet ma- 
terial has drilled a channel through the ambient gas. In this 
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case, the jets of microquasars can propagate along this chan- 
nel without significant deceleration even if the density of the 
jet material is much lower than that of the surrounding gas. 
At the end of their channels the jets decelerate violently due 
to the ram pressure of the ambient gas and the jet material in- 
flates lobes enveloping the jets. This idea w as first propose d 
for the jets of radio-loud active galaxies bv IScheuea (119741) . 
The lobes of these objects give rise to a rich morphology of 
large scale structure detected in radio synchrotron emission 
(e.g. Muxlow & Garrington 1991). Here we suggest that the 
jets of microquasars, and those of GRS 1915+105 in partic- 
ular, give rise to radio lobes emitting synchrotron radiation 
analogous to those obs erved in powerful extra galactic radio 
sources of type FRII ( Fanaroff & Riley 1974). In the case 
of SS 433 radio lobes inflat ed by the jets are indeed observed 
(e.g. lBegelman et alll980h . 

Searches for a similar structure around the first micro- 
quasar identified, GRS 191 5+105, have been unsuccessful 
(Rodriguez & Mirabel 1998, hereafter RM98). No radio syn- 
chrotron structure comparable to the lobes of radio galaxies 
has been found. Only two iras regions with a flat radio spec- 
tra symmetric about GRS 1915+105 and at a distance of sev- 
eral parsec were found. Their flat radio spectra are incom- 
patible with synchrotron emission and they are considerably 
smaller than the expected radio lobes. One possible explana- 
tion for the apparent lack of lobes would be that the jets of 
GRS 1915+105 do not end in strong shocks which acceler- 
ate electrons to relativistic speeds. Only if relativistic elec- 
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trons are fed into the lobes can we expect to detect radio syn- 
chrotron emission. However, it is not clear why the formation 
of shocks should be suppressed at the end of the jets of GRS 
1915+105, since active jets in extragalatic FRII sources al- 
ways end in strong shocks, the so-called radio hotspots. 

In this paper we develop a self-consistent model of the large 
scale structure of the jets of GRS 1915+105. We show that by 
interpreting the two iras regions as impact sites of the jets 
we can successfully apply a fluid dynamical model originally 
developed for extragalactic jets to this source. The model in- 
cludes strong shocks at the ends of the jets, but the resulting 
radio synchrotron emission is too faint to be detected in cur- 
rent radio maps. All other observational results are consistent 
with our interpretation with the exception of one determina- 
tion of the distance of GRS 1915+105. However, we discuss 
this issue in detail in Section |4] 

In Section |2| we review the detailed observat ions of the 
two iras regions by RM98 and Cha tv et all (1200 lL hereafter 
CRM01). We also derive some basic properties of the gas that 
gives rise to the detected radiation. Section |3] describes the 
jet model used and the results obtained from its application 
to the observations. We discuss the implications of the model 
for the distance of GRS 1915+105 and for the velocity of the 
jets in Section 0] Finally, Section |5] summarises our results. 

2. EMISSION FROM IRAS 19124+1106 AND IRAS 19132+1035 

In this Section we discuss the observations obtained in var- 
ious wavebands from the two iras sources. We will show in 
the next Section that all of these observations are consistent 
with an interpretation of the iras sources as sites of excita- 
tion of the ISM heated by the strong shocks driven into this 
material by the GRS 1915+105 jets. 

2. 1 . Emission from hydrogen 

The flat radio spectrum studied in detail by RM98 (see 
also Figure[l} strongly suggests a bremsstrahlung origin from 
ionised atomic hydrogen gas with a temperature of at least 
10 4 K. This view is supported by the detection of the H92a re- 
combination line with a line width of 25 km s" 1 (RM98). As- 
suming that the width of the line is caused by thermal Doppler 
broadening, we would expect a speed of sound of 12.5 km s" 1 
for the emitting hydrogen gas, corresponding to a ga s temp er- 
ature of T = 1.2 x 10 4 K. Note here that in Section ITO we 
will identify the hydrogen gas discussed here with the impact 
region of the jets from GRS 1915+105. The complexity of the 
fluid flow in these regions may well contribute to the broad- 
ening of the H92a line. The temperature of the gas derived 
here is therefore an upper limit. 

The monochromatic luminosity density of ionise d hydro- 
gen d ue to thermal bremsstrahlung is given by (e.g. Longair 

n 2 I hv \ 

L v = 6.8 x 10~ 50 g(v,T)— Vexp -- ergs s" 1 Hz" 1 , 

\ k B T) 

(1) 

where we have assumed that the number density of the elec- 
trons, n, equals that of the protons and no other ions are 
present. We also assume that the density and temperature 
are uniform inside the emitting volume, V. The Gaunt factor 
g (v, T), depends weakly on the temperature and the observing 
frequency, v. 

Given the flux density measurements in the radio, we can 
use equation |T} to estimate the density of the ionised hydro- 
gen gas in the two iras sources. For this we assume that all of 



the emitting gas has the temperature estimated from the width 
of the H92a line. We do not have any information regarding 
the extent of the emission regions along our line of sight or 
their exact three dimensional geometrical structure. Thus it 
is not straightforward to determine the volume of the shocked 
ISM. The continuum radio maps of RM98 show angular ex- 
tents of roughly 35" by 38" for iras 19124+1106 and 40" by 
33" for iras 19132+1035. To account for the uncertainties 
about the source geometry we express their volume as that of 
a sphere with a diameter equivalent to 36.5" at the distance of 
GRS 1915+105 multiplied by an unknown geometrical fac- 
tor /. Given the projected shape of the sources, it is very 
unlikely that / will be much greater than unity. In principle 
the emission regions could be hollow shells seen in projec- 
tion akin to supernova remnants, which would imply / « 1. 
However, in this case we would expect the iras sources to be 
edge-brightened which is not the case. Therefore we assume 
that / ~ 1 . Note also that the volume as well as most other 
quantities derived further down will depend on the distance D 
(measured in kpc) of GRS 1915+105 from us. 

The results for the density of the ionised hydrogen are 
summarised in Table [I] The excellent agreement between 
the density values derived at the three frequencies for each 
emission region supports the interpretation of the emission as 
bremsstrahlung. The densities of the gas in both iras sources 
on either side of GRS 1915+105 are also in remarkably good 
agreement. The table also lists the pressure in this gas for the 
assumption of ideal gas conditions. 

2.2. Molecular line emission 

CRM01 detected sub-mm line emission due to 12 CO, 13 CO, 
H 13 CO + and CS from both iras sources and SiO line emis- 
sion from iras 19132+1035. The existence of molecules in 
both iras sources alone already demonstrates that some of the 
gas must be considerably cooler than the ionised atomic hy- 
drogen. Line widths are only available for SiO (4.4 km s" 1 ) 
for iras 19132+1035. For iras 19124+1106 line widths are 
reported for 12 CO (4.05 km s" 1 ), two velocity components of 
13 CO (3.88 km s" 1 and 3.9 km s" 1 ) and two components of CS 
(3.039 km s" 1 and 1.589 km s" 1 ). If these line widths are due 
to thermal Doppler broadening, then they imply gas tempera- 
tures of a few to several 100s of Kelvin, considerably cooler 
than the atomic hydrogen gas. Both CS and SiO emission are 
considered to indicate high ga s densities comparable to the 
central parts of Hn regions (e.g. Rohlfs & Wilson 2004). The 
low temperatures and high density of the molecular gas prob- 
ably indicate that the iras regions contain a mixture of hydro- 
gen gas with embedded denser clouds of molecular material. 
It is also interesting to note that the shock tracer SiO is only 
detected in iras 19132+1035 but not in iras 19124+1106. 
However, it is not clear from the data presented in CRM01 
whether the molecular material is indeed co- spatial with the 
ionised hydrogen gas. 

2.3. Synchrotron emission from gas 

The radio maps of iras 19132+1035 show a feature with a 
steep radio spectrum to the northwest of the main flat spec- 
trum region (RM98). This feature is elongated with a length 
of 17.4" and a width of 7.9". The spectrum is consistent with 
a power law F v oc y a with a = -0.8. The spectral index 
is typical for synchrotron emission and therefore implies the 
presence in this feature of a relativistic plasma consisting of 
magnetic field and a population of relativistic electrons. We 
can use the standard minimum energy argument to estimate 
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TABLE 1 

Properties of the ionised hydrogen gas observed in thermal bremsstrahlung emission. Radio flux density measurements are from RM98. D is the distance of 

THE EMISSION REGION FROM US MEASURED IN KPC AND / IS THE VOLUME FILLING FACTOR WITH RESPECT TO A SPHERICAL EMISSION REGION (SEE TEXT FOR DETAILS). 



the properties of this relativistic plasma (e.g. lLongairll 19941) . 
Here we neglect any possible contribution to the energy con- 
tent of the emission region by protons. For an emission region 
of monochromatic luminosity L v and volume V the strength of 
the magnetic field is given by 



-G (a) — 
2 V 



2/7 



(2) 



The constant G (a) depends on the observing frequency v, the 
spectral index a and to a lesser extent on the cut-offs of the 
synchrotron spectrum, v m i n and v max . We do not know over 
which range in frequency the synchrotron spectrum extends. 
We therefore arbitrarily set Vmi n = 10 7 Hz and v max = 10 11 Hz. 
Changing either of these parameters does not change our re- 
sults significantly. 

The radio observations of the non-thermal feature imply 
a luminosity density of L v = 6.0 x 10 18 D 2 ergs s" 1 Hz -1 at 
1.5 GHz. We assume a cylindrical geometry for the fea- 
ture with a radius equivalent to 3.9" at the distance of GRS 
1915+105 and a length of 17.4" (sin 6)~\ where 6 is the an- 
gle of the long axis of the cylinder to our line of sight. Using 
equation Q we find the strength of the magnetic field, B^ m . 
For minimum energy conditions the total pressure in this fea- 
ture is then given by 
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(3) 



where k is the ratio of the internal energy stored in any par- 
ticles which do not contribute to the emission of synchrotron 
radiation to the sum of the energy in the magnetic field and in 
the relativistic electrons. All results are summarised in Table 

2.4. Dust 

The two IRAS sources detected by CRM01 (see their Table 
2) show IR emission with a spectral energy distribution of the 
form 5 y ocy- a atv>3x 10 12 Hz (A < lOOyi/m) with a < 3. 
Since the high frequency cutoff in a thermal dust spectrum is 
much steeper than this, we infer a range of dust temperatures. 
The peak of the dust emission must lie at wavelengths longer 
than those observed by CRM01, i e. > 100//m, and therefore 
the bulk of the dust has a temperature of 80K or lower. 

Dust emission is usually modelled using a modified black- 
body spectrum, B(v, T d ), of the form: 

L v = 47TK v B(v,T d )M d . (4) 

The opacity is parametrized as k v oc where ft usually lies in 
the range 1 < J3 < 2. This relationship is normalised such tha t 
at 850/mi (353 GHz), k v = 0.77 cm 2 g" 1 (iDunne et al.l l2000). 
Here we take the value /3 = 2, which g ave the best fit to the ob - 
served properties of nearby galaxies (iDunne & Eales 2001h . 



Frequency (Hz) 

Fig. 1 . — Continuum emission spectrum of the two iras regions from ra- 
dio to near infrared wavelengths. The observational data (symbols) was 
taken from CRM01: Closed symbols: iras 19124+1106, open symbols: 
iras 19132+1035. The lines show our model for the IR emission of iras 
19124+1106 by dust with a range of temperatures (see text for details). The 
dotted lines show the individual contributions to the emission from dust with 
a single temperature. The solid line shows the total emission from all dust 
components. 



and can therefore be used as a good approximation to the av- 
erage properties of the interstellar medium in normal galaxies. 
This value of is consistent with measurements of dif ferent 
environments within our own Galaxy (see e.g. Gordon 1988; 
Sodroski ^eTaliri997t iRistorc elli et a lll998h . 

We model the dust mass distribution as a powerlaw dM = 
C dmt (T/TlT a dT, within the range Tl < T dust < T2. We 
take Tl = 20K, w hich is the average cold dust temperature 
of nearby galaxies (Dunne et al. 2000). The near IR /- and 
K-band observations (at 1.1 and 1.6//m) provide a clear cut- 
off in the emission spectrum, giving a maximum temperature 
T2 = 600K for iras 19124+1106 and T2 = 750K for iras 
19132+1035. The highest temperature components however 
contribute the least to both the total mass and the total energy 
output of the system. In this way, we find good agreement 
with the IR data (see Tableland Figure[T]). 

3. THE JET MODEL 

3.1. Basic setup 

For interpreting the interactions of the large-scale jet struc- 
ture of GRS 1915+10 5 with its environment we will use 
the model developed in Kaiser & Alexander ( 1997, hereafter 
KA97) for the jets of radio galaxies and radio-loud quasars. 



gWmJy V/pc^ jWgg p/10- 10 ergscm- 3 

5 9.9 x 10" 5 D 3 (siney 1 120 (/T 1 sin6>) Z// 4.3 (/T 1 sin6>) 4// (£ + D 

TABLE 2 

Properties of the non-thermal feature next to iras 19132+1035. Radio flux density measurement at 1.5 GHz from RM98. D is the distance of the emission 

REGION FROM US MEASURED FN KPC, 9 IS THE ANGLE OF THE JET LEADING TO IRAS 19132+1035 TO OUR LINE OF SIGHT AND k IS THE RATIO OF THE ENERGY OF NON-RADIATING 

PARTICLES AND THAT OF THE RELATIVISTIC PLASMA (SEE TEXT FOR DETAILS). 



T2/K a L/lO^ergss- 1 M tot /M Q C dust /M Q K~ l 

19124+1106 600 7.2 3.1 D 2 0.38D 2 0.12D 2 
19132+1035 750 7.8 3.5D 2 0.56D 2 0.19D 2 

TABLE 3 

Properties of the dust associated with the iras emission regions. D is the distance of the emission regions from us measured in kpc. 



Non-thermal emission region: 
Hot spot with jet ending in 

strong shock 1 




Thermal Bremsstrahlung 
emission region described 
in Table 1 

Fig. 2. — Sketch of the model setup overlayed on the 20 cm radio contours of 
iras 19132+1035 by RM98. Dashed lines show model features not directly 
observed in the radio. 



The model is based on twin jets emerging from the central 
source in opposite directions. The supersonic jets end in 
strong shock fronts at the location where they impact on the 
ambient gas. The gas coming up the jets inflate two lobes 
which are overpressured with respect to the environment. The 
lobes therefore expand sideways and thus a bow shock is 
driven into the surrounding ISM. The pressure in the lobes 
confines the jets which are propagating through this region. 
By inflating the lobes the jets create their own high pressure 
but low density environment protecting them from disruption 
due to turbulent gas entrainment. This ensures that the jets do 
not widen much on their way through the lobes and thus the 
momentum and energy transported by the jets are delivered to 
the jet shocks without significant losses. 



Here we identify the two iras sources discussed above with 
the impact sites of the two jets of GRS 1915+105. The ba- 
sic setup is sketched in Figure |2| for iras 19132+1035. The 
postulated large scale jet coming from GRS 1915+105 ends 
in a strong shock which gives rise to a so-called radio hotspot. 
This feature is identified with the non-thermal emission re- 
gion at the NW edge of iras 19132+1035. The jet material 
then inflates a lobe which is not detected at currently available 
sensitivity levels. The impact of the jet and the supersonic ex- 
pansion of the lobe drive a bow shock into the surrounding 
ISM. The compressed and shock heated ISM in front of the 
jet is detected as the iras source. We envisage a similar ge- 
ometry for iras 19124+1106. In the following we will refer 
to the jet leading to iras 19124+1 106 as the northern jet while 
the southern jet is pointing towards iras 19132+1035. 

Another prerequisite of the model is that the ejection di- 
rections of the jets remain reasonably constant. The model 
requires that jet ejection s propagate within the lobes created 
by previous jet ejections. Rodriguez & Mirabel ( 1999) report 
a possible precession angle of the small scale jets of GRS 
1915+105 of about 10°. The two iras regions subtend an an- 
gle of about 34" in the direction perpendicular to the jets. If 
this extent is determined by the precession of the jets on large 
scales, causing a shift of the impact region in time, then a 
precession angle of about 1° is implied. In any case, the pre- 
cession angle is small enough for the jets to remain well inside 
the lobes. 

The model assumes a constant energy transport rate, go, for 
each jet. This transport rate is averaged over the lifetime of 
the jets. It is not identical to the energy transport rate esti- 
mated for the jet components observed on small scales. The 
mass density distribution of the ambient gas is modelled by 
a power law, p = p (r/ao)~^, where r is the distance from 
the centre of the distribution and ao is a scale length. Be- 
cause GRS 1915+105 is situated in the ISM of the Galactic 
plane, we will in the following assume that the ambient gas 
has a constant density po and therefore ft = 0. Thus we do 
not have to specify a value for a^. KA97 show that under 
these circumstances the evolution of the large-scale structure 
created by the jet, i.e. the lobe and the bow shock in front of 
it, is self-similar once the jet e xtends beyo nd its characteristic 
length scale given by (see also iFallell 1991) 



P 2 oQo 
Ml 



-1/4 



Go 



1/4 



(5) 



where Mo is the mass transport rate of the jet. The jet material 
travels with a Lorentz factor y y The second equality then 
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follows from the assumption that the energy transported by 
the jet is the kinetic energy of the jet material only. In the 
following we estimate an upper limit for Lo for the jets of 
GRS 1915+105 and show that the measured lengths of the 
two jets exceed Lo. Thus the model of KA97 is applicable to 
the case of GRS 1915+105. 

In Section lTTl we derived an estimate for the density of the 
gas that has been shock heated by the impact of the jets. This 
gas has been compressed by the action of the bow shocks in 
front of the jets. For a strong shock we expect the density 
to be in creased by a factor 4 com pared to its pre- shock value 
po (e.g. Landau & Lifshitz 1987). Taking this into account 
we derive po under the assumption of pure hydrogen gas (see 
Table©. 

The energy transport rate of the jets has so far only been 
estimated during times when jet components are detectable 
milliarcseconds away from the central region in the radio. 
During such periods, the jets are either (i) quasi-continuous 
with a flat radio to IR spectru m, possibly indica ting self- 
absorbed, conical jets (e.g. [Dhawanetal. 2000; but see 
ICollins et alJl2 002). or (ii) discrete, superluminal ejecta emit- 
ting optically thin synchrotron radiation are detected (e.g. 
Mirabel & Rodrigue3lji 



llFenderet al.Hl 999h. During the 
quasi-continuous phases Fender & Pooley (2000) estimate 
energy transport rates between several 10 37 ergss _1 and al- 
most 10 42 ergss _1 , depending on the composition of the jet 
material and its velocity. Similar values are found for the dis- 
crete ejecta (e.g. Fender ^et all 19991) . Note that these estimates 
depend on the unknown distance D of GRS 1915+105 and 
are valid only for limited periods of time when the jets are 
detectable on small scales close to the central source. As dis- 
cussed above, the radio emission from the jets changes sig- 
nificantly over time and this radio variability may well im- 
ply v ariable energy input into the jets (e.g. Klein- Wolt et alJ 
2002). Furthermor e the fact that GRS 1915+ 105 was only 
discovered in 1992 ( Castro-Tirado et al. 1992) in X-rays im- 
plies that the source goes through activity cycles which may 
also apply to the production of jets. Given the substantial ev- 
idence for the correlated behaviour of jets and accretion discs 
in X-ray transients, we would certainly expect this to be the 
case. For the properties of the large-scale structure formed by 
the jets only the energy transport rate of the jets averaged over 
timescales comparable to the fomation timescale of the large- 
scale structure is important. We will argue in Section 13.2.21 
that this time averaged energy transport rate may be as low as 
10 36 ergs s" 1 . 

Finally, to estimate the characteristic length Lo we need the 
Lorentz factor of the jet material. Estimati ng the Loren tz fac- 
tor of jet flows is notoriously difficult (e.g. Fender 2003). Fur- 
thermore, the observed motion of radio 'knots' in the jets may 
simply reflect the motion of shocks along the jets rather than 
the bulk motion of the jet material. Any estimate of y } will 
also depend on the distance to the source and the orientation 
of the jet to our line of sight. Here we use the conservative 
value of yj = 1.1 as a lower limit thereby maximising Lo. 

The maximum for the characteristic jet length allowed by 
the plausible ranges for the relevant source parameters gives 
Lo < 0.013 (/D) 1/4 pc. The angular distance of the two iras 
sources from GRS 1915+105, 26' 10" for iras 19124+1106 
and 26' 24" for iras 19132+1035, imply jet lengths, Lj, of at 
least 7.5 D (sin#) _1 pc, where is the angle of the jet axis to 
our line of sight. Clearly, the jets exceed their characteristic 
lengths by at least two orders of magnitude. Thus the model 



of KA97 is applicable to the jets of GRS 1915+105. 

3.2. Hydro dynamic properties of the jets 

3.2.1. Age of the jet flows 

In Section 12.11 we argued that the temperature of the hy- 
drogen gas heated by the bow shock in front of the jets is 
T = 1.2 x 10 4 K. For a strong shoc k in monatomic gas we 
expect (e.g. Landau & Lifshitz 1987) 



_ = —Mf, 
Ti 16 1 
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where T x is the temperature of the unshocked gas and M- x is 
the Mach number of the shock with respect to this material. 
Substituting for M- x we can solve this equation for the velocity 
of the bow shock which is roughly equal to the advance speed 
of the end of the jet driving this shock, i.e. 
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T = 23kms _1 , 



(7) 



where m p is the proton mass. Clearly the bow shock will be 
highly supersonic with respect to the unshocked ISM, justify- 
ing our assumption of a strong shock. The kinematics of the 
gas in the region affected by the bow shock will be complex 
due to the 3 -dimensional structure of the shock. It is there- 
fore not clear whether the entire line width observed is due 
solely to thermal broadening. Projection of the velocity vec- 
tors of some of the shocked gas and additional turbulent mo- 
tion could contribute to the observed line width implying that 
our estimated temperature of 1, 2 x 10 4 K is an overestimate. 
In this case, we would also overestimate the velocity of the 
shock Lj . In the following we will use the value of 23 km s" 1 , 
but the reader is reminded that this may well be a upper limit 
for the shock velocity rather than its exact value. 
From KA97 we take the result that 



and therefore 
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where t is the age of the jet. c\ is a dimensionless constant 
which depends on the thermodynamical properties of the jet 
material. It also depends on the aspect ratio, R, of the lobe 
inflated by the jet, i.e. R is the length of the lobe divided by 
its width. With the measured lengths given in Tabled we find 
the ages of the jets from equation to be of order 10 6 years 
(see Table©. 

We currently do not understand the connection between the 
processes in the accretion disc and the formation of jets in bi- 
nary systems. However, GRS 1915+105 is cu r rently in a state 
with low accretion rates. Belczynski & Bulik (2002) estimate 
that this stage can last for about 10 7 years. If the production 
of jets is associated with the low accretion state of this sys- 
tem, then our age estimate for the jets is consistent with this 
finding. 

3.2.2. Energy transport rate of the jets and the pressure in the 

lobes 

With these results we can determine the energy transport 
rate, Qq, of the jets from equation (|8|). The constant c\ is 
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Northern j et S outhern j et 

Lj/kms" 1 23 23 

Lj/pc 7.5£> (sinfl) -1 7.8£> (sinfl) -1 

po/10" 22 gem" 3 3.1 (Dfr 1/2 2.0 (Dfy 1/2 

t/10 5 years 1.9 Z) (sinfl)" 1 2.0 D (sinfl) -1 

Qo/lO^ergss- 1 2.8£> 3/2 /~ 1/2 (sinA) -2 2.0£> 3/2 /" 1/2 (sinfl) -2 

p/10- 10 ergscm- 3 0.31 CD/) - 1/2 0.20 (Df)~ 1/2 

B^JfiG 22 CD/) - 1/4 (k + 1)" 1/2 18 (DfT l l A (k + 1) -1/2 

S L5 /mJy beam -1 0.14 (aT 7 )^ (sinfl) -2 (k + l) -7/4 0.08 (Df- 7 f /S (sinfl) -2 (£ + 1) 
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TABLE 4 

Properties of the jets and lobes derived from observations (Lj, Lj and po) AND from the model. D is the distance of the iras regions from us measured est kpc, / 

IS THE VOLUME FILLING FACTOR WITH RESPECT TO A SPHERICAL EMISSION REGION, 6 IS THE ANGLE OF THE JETS TO OUR LFNE OF SIGHT, k IS THE RATIO OF THE ENERGY OF 
NON-RADIATING PARTICLES AND THAT OF THE RELATIVISTIC PLASMA AND R IS THE ASPECT RATIO OF THE LOBES INFLATED BY THE JETS (SEE TEXT FOR DETAILS). 



unknown because the lobes have not been detected and there- 
fore their aspect ratios are unknown. However, for extragalac- 
tic jet sources R > 2 is found with very few exceptions (e.g. 
ILeahv & Williamsl 119841 iLeahv et all 119891) . As a limiting 
case we adopt R ~ 2 and thus c\ - 2 in the following. The 
derived values for g for this assumption are given in Table 
01 Note that go depends strongly on the value of c\, but any 
value R > 2 will result in a lower value of go. 

For the usually assumed distance of GRS 1915+105 of 
D = 1 1 kpc and a viewing angle of roughly = 70° the en- 
ergy transport rate of the jets is comparable to the lower end 
of the range of go estimated at tim es when the small scale 
jets are detectable (see Section l3TTl> . This result is consistent 
with the fact that our energy transport rate is averaged over 
the lifetime of the jets, i.e. over timescales of order 10 5 to 
10 6 years, while the small scale jets used for other estimates 
have much shorter lifetimes. In this context it is interesting to 
note that our value for go is not much lower than the energy 
transport rate inferred for the small scale jets. This strongly 
suggests that the active phases during which GRS 1915+105 
produces jets must be at least comparable in length to the pas- 
sive phases when jet production is switched off. On the other 
hand, episodes of very powerful jet ejections are probably rare 
as they do not seem to dominate the time averaged energy 
budget of the jets. 

In the model of KA97 the pressure inside the lobes is given 
by 



p = 0.0675- 



,10/3 



R 2 



PoQi 



,2\ 1 / 3 



L 4 

j 



(10) 



For R - 2 we find that c x /R ~ 1 and the resulting 
lobe pressures are given in Table |4| For the density of 
the unshocked ISM derived above and an assumed temper- 
ature of this gas of 50 K, we find a pressure of roughly 
1.1 x 10" 12 (Df)~ 111 ergs cm -3 . The lobes will clearly be over- 
pressured with respect to the unperturbed ISM, as required by 
the model. 



3.3. Emission properties of the model 

3.3.1. Radio emission from the lobes 

In the case of extragalatic jet sources a substantial fraction 
of the energy transported by the jets gives rise to a popula- 
tion of relativistic electrons and a magnetic field in the lobes. 
The synchrotron emission from this relativistic plasma is eas- 
ily detected in many such objects in the radio. It is usually 



assumed that the relativistic electrons are accelerated by the 
first order Fermi process at the strong shocks at the ends of 
the jets. Even if they ar e transported by the jets al l the way 
from the central source ( Gopal-Krishna et al. 2001), there is 
no a priori reason to assume that the lobes of microquasar jets 
should not also produce radio synchrotron emission. How- 
ever, searches for such lobes have been unsuccessful for GRS 
1915+105 (RM98). In the following we will derive an upper 
limit for the surface brightness of the expected radio emission 
from the lobes and show that it is well below the sensitivity of 
the radio maps produced to date. 

Assuming that the magnetic field is tangled on scales much 
smaller than the lengths of the jets, we will have three con- 
tributions to the overall pressure in the lobe derived in the 
previous section: The pressure due to the magnetic field, the 
pressure of the relativistic electrons and that of any other par- 
ticles in the lobes which do not emit synchrotron radiation, 
e.g. protons. As in Section l23l the parameter k is the ratio 
of the energy density of the non-radiating particles and the 
sum of the energy densities of the relativistic electrons and 
the magnetic field. Since we want to derive an upper limit for 
the radio surface brightness, we assume that the usual condi- 
tions of minimum energy hold in the lobe. In this case, the 
maximum synchrotron luminosity is obtained from a given 
amount of energy (e.g. lLon gair 1994). We can neglect any 
energy losses of the relativistic electrons due to radiation pro- 
cesses because the age of the source derived above is well 
below the timescale for radiative losses. Solving equation 
for the strength of the magnetic field we find B^ m \ see Table 

El 

We can now solve equation Q for the monochromatic lu- 
minosity L v of the lobes predicted by the model. However, 
because of the lack of direct detections of emission from the 
lobes we have to guess the values of v m i n , v max and a. The cut- 
offs of the spectrum do not greatly influence our result and we 
set them to v m i n = 10 7 Hz and v max = 10 11 Hz as in Section 
12.31 We will argue further down that the non-thermal emis- 
sion region attached to iras 19132+1035 is the strong shock 
of the southern jet. The synchrotron emission from this re- 
gion has a spectral index of a - -0.8. Since in our picture 
the relativistic plasma passing through this shock inflates the 
southern lobe, it is reasonable to use the same spectral index 
to derive L y for the lobes. 

A single telescope beam of angular diameter 6 roughly con- 
tains the flux emitted by a source volume given by 

V = tt(DS/2) 2 1, (11) 
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with 

/ = Lj/(tfsin6>) (12) 

the intersection of maximum length of our line of sight with 
the lobe. Note here that R is the real aspect ratio of the lobe. 
For a viewing angle 6 + 90° this differs from the observed 
aspect ratio by a factor sin 0. Solving equation for L v and 
converting this to flux density, we can now determine the up- 
per limit to the surface brightness of the lobe, S v , as predicted 
by the model. The observational limits are set by the data 
of RM98 with an approximate noise 1-cr rms noise level of 
0.2 mJy per beam with S ~ 4" at v = 1.5 GHz. The limits are 
the same for the other, higher observing frequencies. Given 
the expected power law shape of the spectrum we use this low- 
est observing frequency, i.e. 1.5 GHz, in our estimate because 
it provides the tightest constraint. The results are summarised 
in Table |4| Note that for thinner lobes with R > 2 these sur- 
face brightness limits decrease further. Any deviation from 
the conditions of minimum energy will further reduce the sur- 
face brightness. Unless / « 1 and/or our viewing angle to 
the jets is small, this upper limit is below the detection limit 
of current radio observations. This conclusion is almost inde- 
pendent of the distance of GRS 1915+105. Therefore it is not 
surprising that the synchrotron radio emission of the lobes has 

not been dete cted t o date. 

Recently, iHeinzl <2002l) (see also iHeinzl 120031) has devel- 
oped another approach for estimating the surface brightness 
of the possible radio lobes of microquasars similar to the one 
presented here. Although there are difference in detail s, the 
mode ls are sufficiently similar to warrant a comparison. Heinz 
(2002) concluded that to explain the non-detection of any lobe 
emission from GRS 1915+105 it must be located in an envi- 
ronment with an unusually low density of below 10" 3 cm -3 . 
This apparent discrepancy between the models is explained 
by their assumption of an energy transport rate of 10 46 ergs s" 1 
for the jets. Although appropriate for the jets during major 
outbursts with powerful small scale jets, our model suggests 
that this is an overestimate for the time averaged energy trans- 
port rate. Using ou r value of roughly 10 37 ergs s" 1 we find that 
the predictions of Heinz (2002) agree with ours, namely that 
the lobes should not be detectable in current radio maps. 

3.3.2. Radio emission from the hotspots 

The sites of the strong shocks at the end of the jets of 
extragalactic jets are often referred to as radio hotspots or 
simply hotspots. It is usually assumed that the relativistic 
electrons giving rise to the synchrotron emission of the ra- 
dio lobe s of these objects are acce lerated in these hotspots 
(but see Gopal-Krishna et all 120011 for an alternative view). 
Because of the availability of freshly accelerated relativistic 
electrons and the enhancement of the magnetic field in these 
high pressure regions, the hotspots are considerably brighter 
than the radio lobes. If our model is correct, then the jets of 
GRS 1915+105 should also end in strong shocks and there- 
fore produce radio hotspots which should be easier to detect 
than the lobes. The non-thermal emission region at the north- 
western edge of iras 19132+1035 could be such a hotspot. 
The model we use in previous sections does not make direct 
predictions for the radio emission from hotspots. However, 
there is some circumstantial evidence supporting the identifi- 
cation proposed here. 

The non-thermal feature near iras 19132+1035 is clearly 
elongated in the radio maps. The longer axis points from iras 
19132+1035 back towards the position of GRS 1915+105. 



Both RM98 and CRM01 quote this property as the only ten- 
tative evidence that the two iras regions are connected to GRS 
1915+105 by jets. 

The pressure inside the feature derived in Section 12.31 is 
close to the pressure of the shocked ISM (see Tables [l] and 
13. The model presented here is based on the assumption that 
the pressure at the end of the jets is balanced by the pres- 
sure of the surrounding gas compressed by the bow shock. 
Our findings are consistent with this pressure equilibrium and 
therefore imply that the pressure in the non-thermal feature 
is dominated by the contribution of the relativistic gas, i.e. k 
is small. A significant contribution to the total energy den- 
sity by non-radiating particles would imply that the feature is 
overpressured with respect to the shocked ISM. The approx- 
imate equality of the pressure of the non-thermal feature and 
the shocked ISM also makes it unlikely that the relativistic 
plasma in the feature is far from the conditions of minimum 
energy. 

For an approximately cylindrical lobe we expect the ratio of 
the pressure in the hotspot and that in the lobe to be roughly 
4R 2 (KA97). Identifying the non-thermal feature with the 
hots pot of the southern jet we find from the results of Section 
IT3land l3T2l 

R ~ 2.1 ZT 1/28 (sin<9) 2/7 (k + 1) 1/2 f ll \ (13) 

This aspect ratio is comparable to those found for lobes 
of extragalactic j et sources (e.g. Leahy & Williams! 11984 
Leahv ^etalll 19891) and justifies our assumption of R ~ 2. 

At 1.5 GHz the predicted surface brightness of the non- 
thermal feature is at least a factor 250 greater than that of 
the undetected southern lobe. However, the scale of this dif- 
ference is not surprising as the surface brightness of hotspots 
of extragalactic jet sources can be orders of magnitude greater 
than that of their lobes. 

There is no equivalent non-thermal emission region asso- 
ciated with the northern jet. Given that the small scale jets 
are known to be transient, it is possible that the jets on larger 
scales and therefore the hotspots at their ends are not contin- 
uously fed with relativistic plasma. Any relativistic plasma 
transported by the southern jet into the non-thermal feature 
would pass through the feature along the long axis within a 
time 

W = 100 D (vsinfl)" 1 days, (14) 

where v is the velocity of the plasma in units of the speed 
of light. Clearly the hotspots could therefore be transient 
phenomena as well which disperse on timescales of several 
months to a few years. Subsequently they may be unde- 
tectable for similar timescales until more plasma arrives from 
the central source via the jets. 

In this context it is interesting to note that the shock tracer 
SiO was only detected in iras 19132+1035 which also con- 
tains the non-thermal emission feature (see Section l2T2l . If 
jet material is currently arriving at the end of the jet pointing 
towards iras 19132+1035, then this material should form a 
hotspot, which we identify with the non- thermal feature, and 
it should drive a bow shock into the ISM giving rise to the 
SiO line emission. For iras 19124+1106 no jet material is 
currently arriving at the end of the jet and so no hotspot and 
no SiO lines are detected. 

3.3.3. Dust emission 

In our model the bow shock in front of the jets and lobes 
will heat the ISM and any dust mixed in with it. This 
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process will give rise to at least some of the IR emission 
observed in the two iras regions. One obvious question 
arising from this scenario is as to whether it is realistic to 
find dust cospatial with hot gas (T ~ 10 4 K), and in the harsh 
environment of the bow shocks driven by the jets. However, 
the survival of dust grains has been studied extensively, such 
as in regions of jet-gas interactions in active galaxies (e.g. 



Villar-Ma rtfn. M. and De Young. D. and A lonso-Herrero. A. and 
20011) . and in supernova remnants (e.g. iMouri & Taniguchil 
2000). In the case of slow-moving shocks, such as here, the 
destruction of dust grains is f ound to be a comparatively 
inefficient process De Young (1998), and we would not 
expect them to be destroyed by the jets of GRS 1915+105. 

The observed IR luminosities are very high and the emis- 
sion regions seem to extend well beyond the confines of the 
radio bremsstrahlung emission regions for both iras sources. 
These observations may imply an additional mechanism for 
the heating of the dust on large scales, away from the imme- 
diate impact sites of the jets. Also, the energy transfer from 
hot gas to dust with the addition of radiative losses is a compli- 
cated process and modelling this is beyond the scope of this 
paper. However, the emission peak at 15 /mi coincides with 
the peak of the bremsstrahlung emission in both iras regions 
(see Figure 6 of CRM01) suggesting a close link between the 
jet impact sites and the heating of the dust. 

4. THE DISTANCE TO GRS 1915+105, THE VELOCITY OF THE JETS 
AND THEIR ORIENTATION 

4.1. Distance estimates 

4.1.1. Distance implied by the model 

If the two iras sources are indeed the endpoints of the jets 
of GRS 1915+105, then we would expect GRS 1915+105 
and both iras sources to all lie at a similar distance from us. 
RM98 observed the centre of the H92a recombination line 
to be shifted from its restframe frequency by 57.3 km s" 1 for 
iras 19124+1106 and by 75.7 km s" 1 for iras 19132+1035. 
If these shifts arise from the difference in the Galactic orbital 
velocity of these objects with respect to us, then they imply 
distances of between 6kpc and 7.5 kpc. The observed line 
shifts are essentially confirmed by the millimeter molecular 
lines detected by CRM01, but given the narrower lines and 
therefore smaller uncertainties in these observations, their line 
shift of about 67 km s" 1 for iras 19132+1035 is more likely to 
be correct. 

In our model an additional velocity component is added to 
the gas in the thermal bremsstrahlung emission region by the 
action of the bow shock accelerating the gas in front of the 
jet ends. Before we can determine the distance of the iras re- 
gions from us by using the Galactic rotation curve, we need 
to correct for this effect. Gas at rest with respect to the mean 
rotation of the Galaxy local to the gas will show a line shift 
corresponding to the velocity vlsr, where 'LSR' stands for 
Local Standard of Rest. If this gas is pushed by the jets, then 
the gas in front of the jet at angle 6 to our line of sight point- 
ing preferentially towards us, the approaching jet, will have a 
velocity with respect to us equal to 

Vap = V L SR,ap + U COS 6. (15) 

By symmetry, the gas in front of the jet pointing preferentially 
away from us, the receding jet, will have a velocity equal to 

Vrec = V L SR,rec ~ U COS 6. (16) 

Eliminating the term Lj cos 6 between these two equations 
yields an expression relating the velocities VLSR,ap and 

^LSR,rec • 



Both depend only on the known positions on the sky of the 
two gas clouds, the distance of the gas clouds from us and on 
the adopted rotation curve of the Galaxy. From observations 
of the small scale jets we find that the southern jet associated 
with iras 19132+1035 is the approaching jet and so we set 
v ap = 6 7kms~ 1 and = 5 7.3 km s" 1 . Using the rotation 
curve of Brand & Blitz ( 1993) we then find two possible so- 
AllEriDM-.StockiBln^ttefdc. he distance of GRS 1915+105. Either 
we have VLSR,ap = 61.6 km s" 1 and VLSR,rec = 62.8 km s" 1 im- 
plying Di = 5.4 kpc, or v L sR,a P = 61.4 km s" 1 and v L sR,rec = 
63.0 km s" 1 with D^ - 6.5 kpc. Even after correcting for the 
effects of the bow shock the inferred rotation velocities of the 
gas clouds are comparable to the maximum rotation velocity 
close to the tangent point in the Galaxy at a distance of 6 kpc. 
This explains the small difference between D\ and D^. 

We derive the distance estimates above under the implicit 
assumption that the correction factors in equations JT5l and 
fT6l> are equal in magnitude and opposite in sign. The com- 
plex gas flow in the vicinity of the bow shocks in front of the 
jets makes it impossible to find a definite value for the velocity 
of the bow shock, Lj and thus for the 'real' correction factors 
for v ap and v rec . Therefore there is an uncertainty associated 
with the distance estimates D\ and D^ and in particular with 
the orientation angle 6 ~ 80° implied by equations JT5l and 
fT6l> . However, as shown by RM98 the large measured line 
shifts of at least 57.3 km s" 1 rule out a distance of the iras re- 
gions of greater than 7.4 kpc. It is therefore not surprising that 
the distances estimated from the Galactic rotation curve here 
are considerably smaller than the distance usually assumed 
for GRS 1915+105. 

4.1.2. Other distance estimates 

As mentioned in the previous section, the distance estimates 
to the two iras sources do not allow them to be located far be- 
yond the tangent point at 6 kpc. Attempts have been made to 
measu re the distance to GRS 1915+105 itself. Dhawa net alJ 
(2000) used the Galactic Hi absorption observed during ejec- 
tions of radio emitting small scale jets. They conclude that 
GRS 1 9 1 5 + 1 05 must be at a distance beyond the tangent point 
at 6.1 kpc but closer than 12.2 kpc. By comparison with a 
nearby Hn region they su ggest a distance of greater than 9 kpc. 
Chapius & Corbell (120041) study the optical extinction towards 
GRS 1915+105 for a distance estimate. They derive a very 
similar lower li mit to the distance of 6 kpc, but also show that 
the argument of Dhawan et al. ( 2000) for D > 9 kpc does not 
apply. Both studies agree that the source should be beyond the 
tangent point at 6 kpc. This suggests that if the iras sources 
are in fact associated with GRS 1915+105, then our second 
estimate D2 = 6.5 kpc is likely to be more accurate than D\. 

The observed motion of components in the small scale jets 
also provides a constraint on the distance. As these compo- 
nents cannot move intrinsically faster than the speed of light, 
Fend er etalJ fl999) show that GRS 1915+105 cannot be fur- 
ther away than 11.2 kpc. Unfortunately, this method does not 
provide us with a lower limit for D. 

Finally, two attempts have been made to directly measure 
the motion of GRS 1915+105 in the Galaxy and thus de- 
termin e D using the Galactic rotation curve. Dhawa n"et alJ 
(2000) determined the proper motion of the radio counter- 
part of GRS 1915+105. They found a velocity of 5.8 + 
1.5 mas year -1 which is consistent with the limits derived 
from absorption measurements quoted above, i.e. 6 kpc < 
D < 12 kpc. The determination of the system period by 
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Grei ner et alJ d20Qll) provides an estimate of the line of sight 
velocity of GRS 1915+105. The off-set of the measured radial 
velocity curve determines the systemic velocity with respect 
to the observer. From the Galactic rotation curve they then 
find a distance of 12.1 ± 0.8 kpc. Although this measurement 
seems to rule out that GRS 1915+105 and the two iras re- 
gions are at a common dista nce of 6 kpc , it sho uld be borne in 
mind that the main aim of GreineretaDiOOl) was the deter- 
mination of the orbital period of the binary system. However, 
the orbital period can be obtained independent of the absolute 
calibration of the observed spectra and systematic effects may 
well have altered the determination of the systemic velocity. 
An independent confirmation of the off- set of the radial ve- 
locity curve and thus the systemic velocity would clarify this 
issue. 

The model presented here does not directly depend on the 
distance of GRS 1915+105. In other words, our results re- 
main valid even if the source is located at a distance of 1 1 or 
12 kpc from us. However, in this case the large velocity of the 
two iras sources along our line of sight cannot be explained by 
Galactic rotation. In principle the extra velocity could come 
from a kick the system received during the formation of the 
black hole. This would then mean that the entire structure, i.e. 
GRS 1915+105 and both iras regions, are moving together 
with the kick velocity. There are several problems with this 
explanation. Firstly, we would still expect to see a systemic 
velocity of GRS 1915+105 itself comparable to that of the 
iras sources. Secondly, the kick velocity would need to have 
a component of order 50kms _1 along our line of sight. The 
ends of the jets only advance with a velocity of 23 km s" 1 and 
so the formation of a stable j et channel and the lobe stru cture 
would be impossible. In fact, Belczynski & Bulik (2002) sug- 
gest that GRS 1915+105 did not experience a kick associated 
with the formation of the black hole. 

In summary we emphasize that the model presented here is 
only valid if GRS 1915+105 and both iras regions are located 
at the same distance from us. The iras regions are almost cer- 
tainly close to the tangent point at 6 kpc. Further observations 
to ascertain the distance of GRS 1915+105 are needed and we 
predict that such studies will find the source to be located at 
6.5 kpc from us. 

4.2. Jet velocity and orientation 

As mentioned above, studies of the observed motion of 
components of the small scale jets provide us only with an 
upper limit on the distance of the source. Furthermore, with- 
out a secure distance measurement we cannot determine the 
velocit y of the jet material, vj, separately from the jet orien- 
tation. Fende 7a~rt\ C999) derive an expression for the jet 
orientation to our line of sight depending on D (their equa- 
tion 4). From this expression, it is straightforward to infer the 
velocity of the jet material. Using our distance estimates D\ 
and £>2 from above based on the Galactic rotation curve, we 
find for their radio data 6\ =47° and vji = 0.61 c for D\ and 
62 = 53° and Vj2 = 0.68 c for D^. Both velocities imply only 
mildly relativistic fluid flow in the jets with Lorentz factors 
a round 1 . 3 . 

iFended (120031) points out that for significantly relativistic 
jet speeds we are biased to observe jets at their maximum 
distance from us. If the jets of GRS 1915+105 were highly 
relativistic, this would argue for a distance close to 1 1 kpc for 
GRS 1915+105. However, for the distance required by our 
model the jet flow is only mildly relativistic and so the argu- 
ment does not apply. 



The velocity of the jet and its orientation to ou r line of sight 
inferred for a distance of 1 1 kpc by the method of Fender et al. 
(1999) are vj = 0.98c and = 66°. Any emission line 
produced by the jet material at a rest- frame frequency vo 
would therefore be Doppler- shifted to v app = 0.33vo for the 
approaching jet and to v rec = 0.14vo for the receding jet. 
These results are clearly different from the predictions of our 
model of v app = 1.36vo, v re c = 0.56vo in the case of D\ and 
y app = 1.24vo, y r ec = 0.52vo for£>2. The question over the cor- 
rect distance to GRS 1915+105 could therefore be resolved if 
a pair of emission lines from the small-scale jets could be de- 
tected. In practice the large extinction towards the source may 
make such an observation impractical. 

Note, that the orientation angle estimated here does not 
agree with that implied by equations fT5l and fT6l) . As men- 
tioned in the previous section, the constraints on 6 from the 
Galactic rotation curve are not very strong given the uncer- 
tainties for Lj and the appropriate correction factor for the ro- 
tation velocities of the iras regions. 

5. CONCLUSIONS 

In this paper we have applied the model of KA97 for the 
large scale structure of jets in extragalatic radio sources to 
the microquasar GRS 1915+105. The jets emerging from the 
central source end in strong shocks where they impact on the 
surrounding ISM. The jet material subsequently inflates over- 
pressured lobes which expand at supersonic velocities. The 
ISM surrounding the lobes is compressed by the resulting bow 
shock. This shock is strongest in front of the ends of the jets. 

Although no radio emission from the lobes themselves is 
observed, we identify the two iras regions iras 19124+1106 
and iras 19132+1035 as the jet impact regions. The impact 
regions are almost perfectly symmetric about the position of 
GRS 1915+105 and they lie in the directions that the small 
scale jets of this source point to. The radio emission from 
these impact regions has a flat spectrum which we interpret as 
bremsstrahlung from the shock compressed and heated ISM. 
From the flux density of this emission we derive the density 
of the ISM, while the width of the H92a recombination line 
is used to infer the velocity of the bow shock. Together with 
the observed length of the jets, the density of the shocked gas 
and the shock velocity allow us to derive the age and energy 
transport rate of the jets. The latter is lower than the jet power 
estimated for the large, discrete ejection events observed on 
small scales. However, these events are short-lived and the 
energy delivered to the lobes by the jets is governed by the 
energy transport rate averaged over the lifetime of the jets, of 
order 10 6 years. 

From the pressure predicted by the model we derive an up- 
per limit for the radio surface brightness of the lobes under 
the assumption of minimum energy conditions. We show that 
current radio observations have not gone deep enough to re- 
veal the lobes. We do not confirm claims of an unusually low 
density of the sourc e environment based on this non-detection 
(Heinzl l200i 120031) . Instead, the low radio surface bright- 
ness is caused by the intermittent jet activity that leads to a 
lower average energy transport rate compared to major ejec- 
tion events observed on small scales. 

Further support for our interpretation comes from the ob- 
servation of a non-thermal, radio emitting region near to iras 
19132+1035. This feature is elongated and appears to point 
back towards GRS 1915+105. We identify this feature as the 
strong shock at the jet end or 'radio hotspot' of the southern 
jet of GRS 1915+105. The pressure inside this region derived 
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from its radio emission and projected size is consistent with 
this idea. The absence of a similar feature for the northern jet 
pointing towards iras 19124+1106 can be explained by the 
transient nature of radio hotspots and the intermittency of the 
jet activity. 

The observations of molecular lines and dust emission co- 
spatial with the iras regions is consistent with our interpreta- 
tion. 

The model requires that GRS 1915+105 and the two iras 
regions are located at the same distance from us. Emis- 
sion line shifts place the iras regions at a distance of about 
6.5 kpc. This value is consistent with distance estimates for 
GRS 1915+105 from absorption studies in Hi and the opti- 
cal waveband as well as with the observations of motion in 
the small scale jets. The systemic velocity of GRS 1915+105 
was derived as a byp roduct of the determ ination of the period 
of the binary system (Grein eret alJl2001l) . This measurement 
suggests a distance of 12 kpc to the source, also consistent 
with the results of the other methods, but in conflict with our 
model. 

Our model is a self-consistent interpretation of all cur- 
rently available data on the large scale jet structure of GRS 
1915+105 with the possible exc eption of the distance of this 
source from us as determined by Greiner et al. (2001). How- 
ever, the correct value for the distance is still uncertain and 
will require further observations. One possibility to pin down 
the correct distance would be detection of an emission line 
pair from the approaching and receding jets. The Doppler 



shifts predicted for such a line pair are significantly different 
for different distances of the source. Our model predicts that 
such studies will find a distance of about 6.5 kpc. The model 
further predicts that the radio hotspots of the jets should be lo- 
cated close to the two iras regions. Due to the intermittency of 
the jet activity and their transient nature, follow-up radio ob- 
servations of the iras regions should show significant change 
of the morphology of the southern non-thermal hotspot and 
may very well reveal a similar emission region in the vicin- 
ity of iras 19124+1106. Finally, the radio lobes themselves 
should be detectable in deeper radio observations, possibly at 
lower observing frequencies. 

The model presented here is a first step towards an under- 
standing of the impact of microquasar jets on their environ- 
ment. The energy transported by the jets is found to give rise 
to observable features in the ISM. Further studies of objects 
other than GRS 1915+105 are needed to confirm this picture. 
The large scale structure to which the jets of microquasars 
give rise is comparable to that caused by extragalactic jets. 
Thus the way in which the jets interact with their environ- 
ments is another similarity between the two object classes. 
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